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Abstract 

We explore the phase transitions of the ideal relativistic neutral Bose 
gas confined in a cubic box, without assuming the thermodynamic 
limit nor continuous approximation. While the corresponding non- 
relativistic canonical partition function is essentially a one-variable 
function depending on a particular combination of temperature and 
volume, the relativistic canonical partition function is genuinely a two- 
variable function of them. Based on an exact expression of the canon- 
ical partition function, we performed numerical computations for up 
to 10^ particles. We report that if the number of particles is equal 
to or greater than a critical value, which amounts to 7616, the ideal 
relativistic neutral Bose gas features a spinodal curve with a critical 
point. This enables us to depict the phase diagram of the ideal Bose 
gas. The consequent phase transition is first-order below the critical 
pressure or second-order at the critical pressure. The exponents cor- 
responding to the singularities are 1/2 and 2/3 respectively. We also 
verify the recently observed 'Widoni line' in the supercritical region. 



PACS: 03.75.Hh, 51.30.+i 

Keywords: Critical point. Phase diagram, Ideal Bose gas. Van der Waals. 
* Corresponding electronic address: park@sogang.ac.kr 



1 Introduction 



Spinodal curve, by definition, consists of points in a phase diagram where 
the isothermal volume derivative of pressure vanishes [l|t2j : 

dvP{T,V) = 0. (1) 

The curve corresponds to the border line between thermodynamic stable 
and unstable regions, dvP{T,V) < and dvP{T,V) > respectively [3]. 
Moreover, if it exists, the spinodal curve amounts to a first-order phase 
transition under constant pressure. The temperature derivative at fixed 
pressure acting on an arbitrary physical quantity which is a function of 
temperature and volume, is given by the chain rule of calculus: 

d_ 

dT 

On the spinodal curve, the denominator of the second term vanishes, and 
the temperature derivative of a generic physical quantity along the isobar 
diverges [1] . This provides a possible mechanism for a finite system to mani- 
fest genuine mathematical singularities, without taking the thermodynamic 
limit [S]. 

If we fill a rigid box with water to full capacity and heat the box, the 
temperature will increase but hardly the water evaporates. According to a 
well-known standard argument against the emergence of a singularity from a 
finite system that is based on the analytic property of the canonical partition 
function [3], no discontinuous phase transition should occur. Nevertheless, 
opening the lid will set the pressure as constant or at 1 atm, and the water 
will surely start to boil at 100 degree Celsius. The point is that the usual 
finiteness of physical quantities in a canonical ensemble is for the case of 
keeping the volume fixed. Once we switch to the alternative constraint of 
keeping the pressure constant, first order phase transition featuring genuine 
mathematical singularities may arise from a system possessing finite number 
of physical degrees. The question is then the existence of the spinodal curve. 



d_ 

dT 



'drPiryy 

_dvP{T,V) 



d_ 

dV 



(2) 
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In our previous work we investigated the thermodynamic instabihty of 
the ideal non-relativistic Bose gas confined in a cubic box, based on an exact 
expression of the corresponding canonical partition function. The result was 
that if the number of particles is equal to or greater than a certain critical 
value0 which turns out to be 7616, the ideal Bose gas subject to Dirichlet 
boundary condition indeed reveals thermodynamic instability characterized 
by a pair of spinodal curves and consequently undergoes a first-order phase 
transition under constant pressure. The two spinodal curves are identified 
as the supercooling and the superheating of the ideal Bose gas, 

feBT3,„, = T*pX [7r2ftV(2m)]3/5p2/5 

(3) 

where kg is the Boltzmann constant, m is the mass of the particle, P is an 
arbitrary given pressure, and dimensionless constants of which 

the numerical values depend on the the number of particles, A^, as follows [1]: 



N 




T*p* 


7616 


1.054 3694113 


1.054369 4116 


10^ 


1.052 70 


1.052 77 


10^ 


1.0410 


1.0424 


10^ 


1.034 


1.036 



Between the two temperatures on the isobar, every physical quantity 
zigzags or becomes triple valued, implying the existence of three different 
phases during the discontinuous phase transition. This includes the volume 
as well as the number of particles occupying the ground state. Hence, the 
corresponding first-order phase transition coincides with the Bose-Einstein 

^The number 7616 can be regarded as a characteristic number of a cube. For a sphere 
we get the critical number, 10458. 
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condensation (BEC) in both the coordinate and the momentum spaces. This 
result is comparable to the well known BEC temperature of the ideal non- 
relativistic Bose gas in the thermodynamic limit: 

(4) 

where C is Riemann zeta function. In standard textbooks, e.g. [31EH2], this 
temperature is typically obtained by treating the ground state separately 
and applying the continuous approximation to all the rest of the excited 
states. From the comparison of Eq.(I3]) and Eq.Q, it appears that the con- 
tinuum limit smoothes out some details of the thermodynamic quantum 
structure. 

In the present paper, we generalize our previous work to the relativistic case. 
We investigate into the spinodal curve of ideal relativistic quantum gas of 
neutral bosons confined in a cubic box, subject to Dirichlet boundary con- 
dition. While the non-relativistic canonical partition function is essentially 
a one-variable function depending on a particular combination of temper- 
ature and volume, the relativistic canonical partition function is genuinely 
a two-variable function of them. We show that if the number of particles 
is equal to or greater than the critical number 7616, the ideal relativistic 
neutral Bose gas features a spinodal curve with a critical point where the 
supercooling and the superheating lines converge. Further, we verify the re- 
cently observed 'Widom line' P2] emanating from the critical point as the 
crossover between liquid-like and gas-like behaviour in the supercritical re- 
gion [TT] (see also [T2] and [13]). The critical number coincides with that of 
the non-relativistic system and hence, the present paper confirms our pre- 
vious numerical results 

The main motivation for us to consider the relativistic ideal Bose gas in the 
present paper as a generalization of the non-relativistic system we studied 
previously is, by letting the corresponding canonical partition function be a 
two- variable function, to observe a more intriguing phase diagram such as a 
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critical point and the Widom line. After all, quantum mechanics and rela- 
tivity are two cornerstones of modern physics. For pioneering earlier works 
on the ideal relativistic Bose gas via the continuous approximation we refer 
to Refs. |14fll7j . For recent discussion see [18] and references therein. In 
the present paper, we restrict ourselves to neutral bosons. The inclusion of 
anti-bosons, as in |19H21j . will be treated separately elsewhere |22j . 

The rest of the paper is organized as follows: 

In section [21 we present model independent analysis on canonical ensemble. 
We first review three different expressions for the canonical partition func- 
tion of generic non-interacting identical bosonic particles. We then, simply 
by assuming the existence of a spinodal curve, show that under constant 
pressure first-order phase transitions feature singularities of exponent 1/2 
while second-order phase transitions have the critical exponent 2/3. 
Section [3] is exclusively devoted to the ideal relativistic neutral Bose gas con- 
fined in a cubic box. After analysing non-relativistic and ultra-relativistic 
limits, we present our numerical results. When the number of particles is 
equal to or greater than 7616, the ideal relativistic Bose gas reveals a spin- 
odal curve with a critical point. We draw the corresponding phase diagram. 
The final section [J] summarizes our results with comments. 

2 General analysis 

2.1 The canonical partition function: review 

When a single particle system is completely solvable, each quantum state is 
uniquely specified by a set of good quantum numbers, which we simply de- 
note here by a vector notation, n. With the corresponding energy eigenvalue 
Efi and /3 = l/{kBT), we define for each positive integer a, 

Aa := En e-^^^rt , (5) 
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where the sum is over all the quantum states. In particular, when a = 1, Xa 
coincides with the canonical partition function of the single particle. 



For an A^-body system composed of the non-interacting identical bosonic 
particles above, we may write the corresponding canonical particle function 
in three different manners: 



• As first obtained by Matsubara [23j and Feynman 

Zn = En.. Ua=l {Kr^/{mal oJ^'^) , (6) 

where the sum is over all the partitions of N , given by non-negative 
integers rua with a = 1, 2, • • • , iV satisfying N — E^=i ^ ™a ■ 

• By a recurrence relation first derived by Landsberg |25| . 

^JV- (Er=lAfe^Ar-fc)/iV, (7) 

where we set Zq — 1. 

• From Ref. g], 

Zn = det(r!Ar) (Z^)^ /Nl (8) 

where at is an almost triangularized NxN matrix the entries of which are 
defined by 

Aa-6+i/Ai for h <a 
VLN\a,h\:=\ -a/\i for h = a + l 

otherwise . 

hi particular, every diagonal entry is unity. 

The first expression ^ implies that, compared to ideal Boltzmann gas, 
ideal Bose gas has higher probability for the particles to occupy the same 
quantum state [4j. The last expression is useful for us to see when the 
canonical partition function reduces to the conventional approximation |26] : 

{Zif/N\ (9) 
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which would be only valid if all the particles occupied distinct states, as 
in high temperature limit. Indeed, in this limit we have Ai — >■ oo. Hence 
det(f]7v) — >• 1 and the reduction Q holds. 



According to the Hardy-Ramanujan's estimation, the number of possible 
partitions grows exponentially like and this would make 

any numerical computation based on the expression ([6]) practically hard for 
large N. Among the three above, the recurrence relation d?]) provides the 
most efficient scheme of N'^ order computation. 

2.2 Critical and noncritical exponents: universal results 

When the energy eigenvalue of the single particle system Ef^ ^ depends 
on volume, the canonical partition function depends on temperature and 
volume, so does the pressure: 

P(r, V) = k^Tdv In Zn{T, V) . (10) 

We consider, at least locally|§ inverting P{T, V) to express the tempera- 
ture as a function of P and V like T{P,V). Plugging this expression into 
dyPiT, V), let us define a quantity <!' as a function of P and V: 

HP,V) := - dvPir,V)\^,^^^p^y^ . (11) 

As we change the volume keeping the pressure constant, ^ indicates whether 
the thermodynamic instability develops or not on the isobar. With the minus 
sign in front, now negative ^ corresponds to the thermodynamic instability. 
In a similar fashion to ([2]) we write 

+ ^dpT{P,V)-l^ (12) 



d_ 

dV 



dV 



from which it is straightforward to obtain 

dy^P, V) = -d^P{T, V)-^ dpT{P, V) drdyPiT, V) , (13) 



^Since the vanishing of dvT{P,V) imphes the divergence of 9tV^(-P, T), inverting 
P{T, V) to express the volume as a function of P and T is not conceivable close to the 
spinodal curve, as we see shortly from the relation (|14p . 
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and relate the first and second volume derivatives of T(P, V) to ^{P, V) as 



dvTiP,V) = ^ dpT{P,V) , 

(14) 

d^T{P, V) = dv^{P, V) dpT{P, V) + <^ dpdvT{P, V) . 
Further, on the spinodal curve = 0) - if it exists - we have 

dT _ 9lP{T,V) , . 

dV - dTdvP(T,V) ■ y^'^) 

Now critical point can be defined as follows. At the critical pressure 
P = Pc , the system is generically stable i.e. $ > 0, except at only one 
point that is the critical point. This implies that the critical pressure line 
is tangent to the spinodal curve, such that at the critical point, $ and dy^ 
vanish: 

HPc,Vc) = o, dvHPc,Vc) = o. (16) 

Moreover, from (fTTI) . (fT3]l . (fT^ . these two conditions are equivalent to 

dvP{n, Vc) = 0, d^Pin, Vc) = 0, (17) 

and also to 

dvT{P„ Vc) = 0, dlT{P,, V,) = Q. (18) 

Eq. (|17p corresponds to the usual definition of the critical point as an inflec- 
tion point in the critical isotherm on a (P, V) plane [3]. With (jl5p . Eq. (|17p 
implies that the critical point is an extremal point on the spinodal curve. 
On the other hand, Eq. (jl8p implies that the expansion of T(P, V) near the 
critical point on the critical isobar starts from the cubic order vaV — Vc (for 
related earlier works see [27] and references therein): 

T(Pc, V)-T, = \{V - Vcfd^TiP,, Vc) + higher orders . (19) 

Clearly this leads to the following critical exponent: 

V/Vc-l ~ \T/T^-lf, /3 = l/3. (20) 
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Similarly, from (jl3p and (|17p , the expansion of ^ on the critical isobar starts 
from the quadratic order in V — Vc such that 

^{Pc,v) ~ {v/Vc-if ~ (21) 

Thus, from ([2]), any physical quantity given by the temperature derivative 
along the critical isobar diverges with the universal exponent 2/3. This in- 
cludes the critical exponent of the specific heat per particle under constant 
pressure: 

Cp ~ |r/r, -Ip", a = 2/3. (22) 

On the other hand, since the specific heat per particle at constant volume 
is finite for any finite system, the corresponding critical exponent is trivial: 

Cy = finite i.e. q = 0. (23) 

Note that our conclusion does not exclude the possibility for Cy to develop 
a peak. We merely point out that the peak will be smoothed out if it is 
sufficiently zoomed in. Further from ()2ip . the inverse of $ gives the critical 
exponent of the isothermal compressibility on the critical isobar: 

Kr:=-V-^dpViPc,T) ^ \T/T,-1\~\ 7 = 2/3. (24) 

Finally, from (|17p . we also obtain an isothermal critical exponent at the 
critical temperature, T = Tc'. 

P/Pc-l ~ \V/Vc-l\^, 5 = 2,. (25) 



On a generic isobar below the critical pressure, P < Pc, the thermody- 
namic instability given by ^{P, V) < appears naturally over an interval in 
volume and hence in temperature as well. The two ends then correspond to 
a supercooling point and a superheating point. At these points $ vanishes 
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and between them <I> is negative, as depicted in FIG. [T] from the Van der 
Waals equation of state as an example|f 



[(P/Pe) + m/Vf] [{V/V.) 



(26) 



\ = W/Tc). 

It is crucial to note that the volume assumes triple values between the super- 
cooling and the superheating temperatures at constant pressure less than Pc ■ 




Figure 1: Van der Waals equation of state 

The dashed lines are constant pressure lines of three different values, and the 
thick solid line is the spinodal curve satisfying <I> = —dvP(T,V) = 0. When 
P < Pc the constant pressure line crosses the spinodal curve twice, at the su- 
percooling point (K,T*) and at the superheating point (K*,P**). Between the 
two temperatures the volume is triple valued. If P = Pc the constant pressure 
line comes in contact with the spinodal curve only once at the critical point 
{Vc,Tc). Otherwise, i.e. P > Pc, the constant pressure line does not undergo 
thermodynamic instability and hence no first-order phase transition arises. On 
the spinodal curve the critical point has the highest temperature. 

''For a modern exposition of Van der Waaals forces, see e.g. [28) . 
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In the derivation of the critical exponents above we tacitly assumed 
dY^{Pc,Vc) / 0. In fact, the equivalence among (fT6|) . (fTT]) and (fTH]) gener- 
alizes up to an arbitrary order, say n : For all /c = 1, 2, • • • , n, the vanishing of 
d!^P{T, V) is equivalent to that of 5^T(P, V) as well as that of V). 
Hence, if the first non-vanishing multiple volume derivative of $ occurs at 
the order n as d^^{Pc,Vc) / 0, the exponents in ([20]), ([22]), dM]), (|25]) . 
assume the general values: 

Clearly these satisfy the following "scaling laws" [3ll29j: 
Rushbrooke : + 2/3-^7 = 2, 

(28) 

Widom : 7 = /3(,5 - 1) . 

Especially, the case of n = 1 corresponds to the first-order phase transition 
at noncritical, supercooling and superheating points. Even at these points 
our analysis shows the existence of the universal exponent: 

a = /3 = 7 = 5-1 = 1, (29) 

which agrees with the mean field theory result [30] (see also ^U). Around 
the supercooling or superheating point in FIG.dl t^ie constant pressure lines 
can be approximated by a parabola straight up or down respectively. The 
possibilities of n > 3 shall correspond to exceptional critical points. 

In the next section we show that ideal relativistic quantum gas of neutral 
bosons indeed features a spinodal curve similar to FIG. [TJ 
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3 Ideal relativistic neutral Bose gas 



3.1 Algebraic analysis 

Hereafter, exclusively for ideal relativistic neutral Bose gas, we focus on N 
particles with mass m, confined in a box of dimension d and length L = V^^'^. 
Hard, impenetrable walls impose Dirichlet boundary condition^ 
With positive integer valued good quantum numbers: 

n = (ni,n2, • • • ,nd) , (30) 

the spatial momentum is quantized as p = nhfi/L, such that the single 
particle Boltzmann factor in ^ assumes the form: 

=exp(-r-Vl + n-nV-2/'^ ) . (31) 

Here we have introduced, for simplicity of notation and the forthcoming 
analysis, dimensionless temperature: 

T := k^T/{mc^) , (32) 

as well as dimensionless volume: 

V :=[mc/{^h)Yv . (33) 

The canonical partition function is then two-variable function of them, 
Zn{T,V). We further define dimensionless pressure by: 

V := [(^n)V(m'^+ic'^+2)] P = Tdv In Zn{T, V) , (34) 

and a dimensionless indicator of the thermodynamic instability, as for $ in 

cP := -{l/N)V^dllnZNiT,V) . (35) 

*We recall that nevertheless enforcing periodic or Neumann boundary condition leads 
to thermodynamic instability at low temperature near to absolute zero for arbitrary N [4]. 
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Other physical quantities we are particularly interested in the present paper 
are specific heats per particle at constant volume or under constant pressure: 

Cv = (k^/N) [{rdr)^ + rdr] in ZN{r,v) , 

(36) 

Cp = Cv- k^T[dTV{T, V)]y[NdvV{T, V)] . 



• Non-relativistic limit. The non-relativistic limit corresponds to the limit of 
low temperature and large volume with TV'^^'^ held fixed. For large volume, 
the Boltzmann factor ()3ip reduces to that of the non-relativistic particle, up 
to the shift of the energy by mc^: 

Physical quantities like cj), Cy, Cp become one- variable functions depending 
on the sing le variable T'V^/'^, or alternatively TP~^/('^+^) [4J. The constant 
energy shift is irrelevant to them. Furthermore, at extremely low tempera- 
ture only the ground state energy matters and the canonical partition func- 
tion reduces to 

Zn e-l^^^o = exp(- NT-Wl + dV-^/'^ ) . (38) 
Hence near to absolute zero temperature we obtain 

N/V = V^+^Z'^ ^d + V^ , (39) 

and also 

(j) = oo, Cv = 0, Cp = 0. (40) 

For large volume, Eq. (l39p further leads to V = (A^/J')'^/('^+2) and this agrees 
with Ref . [1] . Apparently the volume assumes a finite value at absolute zero 
temperature, essentially due to the Heisenberg uncertainty principle. It is 
worthwhile to recall 

PV = NksT, (/.= !, Cv/ks = d/2, Cp/k^ = l + d/2, (41) 



12 



in the large TV^^'^ limit of the non-relativistic case |4j. Note that the first 
relation is equivalent to W = NT- 

• Ultra-relativistic limit. The ultra-relativistic limit corresponds to the limit 
of high temperature and small volume with TV^^'^ held fixed. At sufficiently 
high temperature, highly excited states dominate and from ([9]) the canonical 
partition function reduces to (Zi)'^ /N\, where the Boltzmann factor also 
reduces as 

g-/3i?« — ^ ^-Mp\ = exp[-|n|/(rV^/'^)j . (42) 

Consequently, the canonical partition function and the physical quantities 
like (j), Cy, Cp become one- variable functions depending on the single vari- 
able TV^/'^, or alternatively from ([Ml) and (jl2]) on 7"P-V(rf+i). 

Moreover, for large TV^^'^ we may safely assume the continuous approx- 
imation to obtain |32j, 

z^(r,v)oc(vrT, (43) 

which implies in the large TV^^'^ limit, cf. (j4ip : 

PV = NksT, (/) = 1, Cv/ks = d, Cp/ks = l + d. (44) 



3.2 Numerical results 

Here we present our numerical results on the ideal relativistic neutral Bose 
gas confined in a cubic box, i.e. d = 3, at generic temperature and volume, 
performed by a supercomputer (SUN B6048). Our analysis is based on the 
recurrence relation ([7]) and the previous work [^j^ . 



^The numerical computation inevitably requires a momentum cutoff as for n-n, which 
generically deforms the result at high temperature. We choose mutually different, suf- 
ficiently large values of the cutoff, such that we maintain at least for some interval of 
temperature the high temperature behaviour (|44p . and report only the cutoff independent 
results {cf. Refs. [331E1])- 
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• Emergence of a spinodal curve (FIG. [2]). As grows, (p develops a valley 
of local minima, which eventually assumes negative values if > 7616. 




Figure 2: Constant pressure lines and spinodal curve 
Dashed, thin solid or thick solid lines denote respectively the constant pressure, 
constant (p oi' spinodal curve. Near to the origin of the {T,V/N) plane cp 
diverges and in the opposite infinite limit cf) converges to unity. When A^ = 1 
(FIG. [2^ with Vl < Vm < Vh), (p is monotonically decreasing from oo to 1 on 
arbitrary isobars. As A^ increases (p develops a valley whose height is less than 
unity. Moreover if A^ > 7616 the valley assumes negative values and a spinodal 
curve emerges. FIG.[2tl magnifies the tip of the spinodal curve for A^ = 7616 to 
manifest a critical point, {Tc,Vc/N) = (3.6651475 x 10-^ 9.510401 x 10^). 
In FIG.Eb, Vc denotes merely the numerical value, Vc = 2.0151967 x 10"^\ 
which amounts to the critical pressure in the system with one more particle, 
A^ = 7616. 
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• Specific heat per particle under constant pressure (FIGs. [3llll[5]). During 
the first-order phase transition under constant pressure which is less than 
the critical value, the volume and hence every physical quantity are triple 
valued between the supercooling and the superheating temperatures. We 
focus here on Cp, the specific heat per particle under constant pressure. 
For the explicit behavior of other quantities we refer to our earlier work [5]. 



a. N=l h. N = 7615 




Figure 3: 31) image of Cp on {T,V/N) plane for = 1 and N = 7615 
The plane decomposes into three parts: phase I with Cp ~ 0, phase II with 
Cp ~ 2.5A;b and phase III with Cp ~ Ak-^. When = 1 the transitions are 
monotonic and smooth. As N grows, at the borders between I and II as well 
as between I and III, a range of peaks emerges which will eventually diverge 
for N > 7616. FIG. [3b has been cut at the height of Cp = Bk^, and the actual 
peak rises up to Cp ~ 5.336 84 x lO^k^. 
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c. W = 7616 
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Figure 4: Cp for various N and V 

For = 1, Cp is monotonically increasing from zero to on any isobar. 
Under sufficiently low pressure, it assumes the intermediate value of 2.5/cb as in 
Eq. (l4T]l . As N increases, Cp develops a peak on each isobar. Especially when 
N > 7616 and V < Vc, it diverges both to the plus and minus infinities at the 
supercooling point T = % as well as at the superheating point T = %*. At 
the critical point T = % = 3.6651475 x 10^^, it diverges only positively. For 
V > Vc, the specific heat features a single finite peak which can be identified 
as the Widom line. 
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a. supercooling 

IniCp/ks) Ini-Cp/kB) 




b. superheating 

IniCp/ks) \n{-Cp/kB) 




Figure 5: Critical and noncritical exponents 

Our numerical data confirm the singular behavior of the specific heat Cp an- 
ticipated in Sec. 12. 21 above the supercooling point, Cp ~ ±(7" — 7^)^^^^ 
(FIG. [Si); below the superheating point, Cp ~ ±(r** - T)~^^^ (FIG.Eb); 
and around the critical point, Cp ~ IT — Tc\~'^^^ (FIG.Eb)- The numer- 
ical data are for iV = 7616, V = 5.00017005640 x 10"^^ (FIG.[3i,b) or 
Vc = 2.0151967 X 10"^^ (FIG.Eb)- The error bars originate from the nu- 
merical uncertainty beyond these digits. The straight lines correspond to the 

theoretical slopes, -1/2 (FIG.[3i,b) or -2/3 (FIG.Efc). 
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• Phase diagram of the ideal relativistic neutral Bose gas for N = 10^. 
(FIG.ED. 




Figure 6: Phase diagram for N = 10^ on {V,T) and {T,V/N) planes 
Thick solid line, dashed line and dotted line correspond to the spinodal curve 
(including its inside), the Widom line (a range of finite peaks in Cp which 
also coincides with the valley of 0) and the T ~ 1 line respectively. The 
lines divide the phase diagram into three parts: phase I of Cp ~ 0, phase 
II of Cp ~ 2.5A;b and phase III of Cp ~ Ak^. On the top of the solid line, 
there exists a critical point. As N increases, the critical point moves toward 
more ultra-relativistic, higher temperature region along the Widom line keeping 
Tc{Vc/Ny/^ and TcVc^^^ constant. 
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Our numerical results of the critical point are listed below for some 
selected N: 



N 


Tc 


Va/N 




7616 


3.6651475x10-5 


9.510 401 xlO^ 


2.0151967x10-" 


10^ 


0.149 54 


3.133 


0.02702 


10^ 


1.738 


0.025 2 


53.1 



which gives 



N 


Tc{Vc/N)y^ 








7616 


0.354 4519 


0.695 6079 


3.604 329 X 10"^ 


1.729 865 X 10-2 


10^ 


0.320 2 


0.634 


0.2188 


0.368 8 


10^ 


0.149 


0.355 


0.510 


0.644 



Comparison with the non-relativistic results in Ref. [2] indicates that the 
critical point of = 7616 is in the non-relativistic region while those of 
N = 10'' and N = 10^ are not. 
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4 Summary and comments 

In summary, as seen from FIG. [3] and FIG. [H the ideal relativistic neutral 
Bose gas divides the {T,V/N) phase diagram into three parts. 

• Phase I: condensate with Cp ~ (I40p . 

{ (r, v/N) I T < (iv/v)^/^ & r < (iv/v)2/3} . 

• Phase II: non-relativistic gas with Cp ~ 2.5A;b (|4ip . 

{{T,v/N) I r < 1 & r> (iv/v)^/^}. 

• Phase III: ultra-relativistic gas with Cp ~ 4A;b (|44p . 

{ (r,V/iV) I T>{N/Vf/'^&LT>l]. 

Equivalently the ("P, 7~) phase diagram splits into three parts. 

• Phase I: { (P, T) | T < P^/s ^ 7- < pi/4|_ 

• Phase II: [ {V ,T) \ T >V'^'^ k T <l]. 

• Phase III: [{V,T)\T> V^'^ & T > l}. 

Here < and > mean rough inequalities up to constants of order unity. 

When = 1, the transitions are all smooth and monotonic along any 
isobar. As N increases, a valley in </> and hence a range of peaks in Cp 
develops along the boundary of the condensate phase I. Especially, when 
N = 7616, the valley of (j) assumes negative values and a spinodal curve with 
a critical point emerges along the boundary of the phase I starting from 
the absolute zero temperature. Beyond the critical point it is the Widom 
line that divides the phase I and the phase III. As N further increases, 
the critical point moves along the Widom line toward the ultra-relativistic 
higher temperature region. For the transition between the non-relativistic 
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and ultra-relativistic gas phases II and III there is no latent heat involved. 

The spinodal curve sharply defines the phase diagram. The consequent 
phase transition is first-order below the critical pressure or second-order at 
the pressure. The exponents corresponding to the singularities are 1/2 and 
2/3 respectively. Presence of both the supercooling and the superheating 
characterizes the first-order phase transition. 

The resulting equation of state from the spinodal curve of the relativistic 
ideal gas resembles the Van der Waals equation of state which is derived by 
assuming both repulsive potential at short distance and attractive potential 
at long distance. It is well known that the effective statistical interaction of 
ideal Bose gas is attractive |3] . Our result seems to suggest a more intriguing 
structure. The fact that the volume is finite at absolute zero temperature 
seems also to indicate a repulsive interaction at short distance related to the 
Heisenberg uncertainty principle. 

An interesting open question is whether the critical point converges or 
not in the phase diagram as — t- oo. If not, one may also wonder whether 
there exists another critical number in N for the emergence of a spinodal 
curve from the ultra-relativistic canonical partition function ()42p . We have 
verified for up to A = 10^ that the ultra-relativistic canonical partition 
function does not feature any spinodal curve, though develops a local 
minimum. 

Although in this work we have focused on the relativistic generaliza- 
tion of the non-relativistic ideal Bose gas and have obtained a nontrivial 
phase diagram with a critical point and Widom line, it is natural to expect 
that other generalizations which involve extra dimensionful parameters, e.g. 
trapping potentials, may also lead to qualitatively similar or even more in- 
triguing phase diagrams. 
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